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The possibility that dark matter, whose existence is inferred from the study of the galactic rotation 
curves, and from the mass deficit in galaxy clusters, can be in a form of a Bose-Einstein Condensate, 
has been extensively investigated lately. In the present work, we consider a detailed analysis of the 
astrophysical properties of the Bose-Einstein Condensate dark matter halos that could provide clear 
observational signatures that help discriminate between different dark matter models. In the Bose- 
Einstein condensation model dark matter can be described as a non-relativistic, gravitationally 
confined Newtonian gas, whose density and pressure are related by a polytropic equation of state 
with index n = 1. The mass and gravitational properties of the condensate halos are obtained in a 
systematic form, including the mean logarithmic slopes of the density and of the tangential velocity. 

The lensing properties of the condensate dark matter are investigated in detail. In particular, a 
general analytical formula for the surface density, an important quantity that defines the lensing 
properties of a dark matter halos, is obtained in the form of series expansions. This enables arbitrary- 
precision calculations of the surface mass density, deflection angle, deflection potential, and of the 
magnification factor, thus giving the possibility of the comparison of the predicted lensing properties 
of the condensate dark matter halos with observations. The stability properties of the condensate 
halos are also investigated by using the scalar and the tensor virial theorems, respectively, and 
the virial perturbation equation for condensate dark matter halos is derived. As an application of 
the scalar virial theorem we consider the problem of the stability of a slowly rotating and slightly 
disturbed galactic dark matter halo. For such a halo the oscillations frequencies and the stability 
conditions are obtained in the linear approximation. 

PACS numbers: 98.62.Gq; 98.62.Sb; 95.35.+d; 03.75.Nt 


I. INTRODUCTION 


The polytropic gas model, in which the pressure P 
has a simple power law dependence of the matter energy 
density p , given by P = Kp ^, plays an important role 
in the modelling of self-gravitating systems of compact 
astrophysical objects, such as white dwarfs and neutron 
stars, as well as in the study of the equation of state 
of nuclear matter. Usually the polytropic equation of 
state is represented as P = Kp 1+1 / n , where n is called 
the polytropic index. The Poisson and the hydrostatic 
equilibrium equations for a self-gravitating, spherically 
symmetric polytropic fluid can be combined to provide 
the basic equation describing the equilibrium properties 
of the system as 


1 d ( ,, dO \ 

edc v c dc) 


+ 9 n = 0, 


(1) 


where the dimensionless quantities 6 and £ are defined 
as p = p c 0 n and ( = \J\-nGp\l (n + 1) P c r , with p c and 
P c the central energy density and pressure, respectively. 
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Eq. (jTj) is known as the Lane-Emden equation, and its 
physical, astrophysical and mathematical properties have 
been intensively studied by using both analytical and nu¬ 
merical methods. It admits only three exact solutions, 
corresponding to n = 0, n = 1 and n = 3, respectively 

B- 


From these three solutions, the unusual solution corre¬ 
sponding to n = 1, having the form 9 = sin<C/((, does not 
seem at first sight to have many realistic physical prop¬ 
erties, and consequently it has been less investigated, as 
compared to the other solutions. However, somehow un¬ 
expectedly, the n = 1 polytropic equation of state ap¬ 
pears in a fundamental field of research apparently un¬ 
related to astrophysics, namely, the theory of the con¬ 
densed Bose-Einstein cold atomic gases. A basic result in 
quantum statistical mechanics is that at very low temper¬ 
atures, in a dilute Bosonic system (Bose gas) all particles 
condense to the same quantum ground state, forming a 
so called Bose-Einstein condensate (BEC). From a physi¬ 
cal point of view a BEC is characterized by a sharp peak 
in both coordinate and momentum space distribution. 
The BEC phase transition takes place when the particles 
are correlated with each other quantum, mechanically. 
This happens when their wavelengths overlap, that is, 
the thermal de Broglie wavelength A t becomes greater 
than the mean inter-particles distance l. Hence the BEC 
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transition occurs at a temperature 


T < 


27 tH 2 
rnks 


,2/3 


( 2 ) 


where m is the mass of the particle in the condensate, n 
is the number density, and ks is Boltzmann’s constant 
1-i- A coherent quantum condensate state develops in 
two physical situations: a) when the particle density is 
high enough, or b) the system temperature is sufficiently 
low. The Bose-Einstein Condensation has been observed 
in laboratory experiments [M2- A large number of 
quantum degenerate gases have been created by a com¬ 
bination of laser, magnetic and evaporative cooling tech¬ 
niques. The experimental investigations have opened sev¬ 
eral new directions of research, involving atomic, statis¬ 
tical and condensed matter physics 0-01- 

The purpose of this paper will be applications of the 
BEC to the dark matter problem. Indeed, the existence 
of dark matter is primarily supported by the rotation 
curve of galaxies [13| 15], and also by the difference be¬ 
tween the masses derived from luminosity and from the 
virial theorem in clusters of galaxies [16] . Furthermore, 
gravitational lensing data also indicates that light deflec¬ 
tion is more accentuated than would be to visible matter 
alone E3- Thus, by introducing a new matter compo¬ 
nent, denoted dark matter, it is possible to explain all of 
these observational phenomena. However, the nature of 
dark matter, and its constituent particle(s), still remain 
essentially unknown, despite several decades of intensive 
research. For reviews on the dark matter problem see 
00 . 

In this context, the possibility that dark matter could 
be in the form of a Bose-Einstein condensate was ana¬ 
lyzed in detail in [19] (for earlier work in this field see HJ- 
|23|). In particular, the equation of state of the conden¬ 
sate dark matter was obtained as a polytropic equation of 
state of index n = 1, and the corresponding density pro¬ 
file was discussed. The theoretical rotation curves were 
fitted with the observational data from a sample of Low 
Surface Brightness (LSB) galaxies, and it was pointed 
out that the lensing properties may discriminate between 
the standard pressureless and the condensate dark mat¬ 
ter models. The physical and astrophysical properties of 
Bose-Einstein Condensate dark matter halos were inves¬ 
tigated in [24[ . 

Recently, the rotation curves in the Bose-Einstein con¬ 
densate dark matter model were analyzed in [20], and a 
good agreement with observations was found. A com¬ 
parison of the predicted rotation curves with the obser¬ 
vational data for eight dwarf galaxies was performed in 
[20, and it was shown that the presence of the conden¬ 
sate dark matter can solve the core/cusp problem faced 
by the standard ACDM (ACold Dark Matter) model. 
The effects of the angular momentum and on the vor¬ 
tices on the equilibrium of self-gravitating, rotating BEC 
haloes which satisfy the Gross-Pitaevskii-Poisson equa¬ 
tions were considered in [271430 . In a Bose-Einstein con¬ 
densate vortices form as long as the self-interaction is 


strong enough, which is not the case for axionic dark 
matter. The galactic masses in the framework of the con¬ 
densate model were discussed in [30, HU, and the effects 
of the finite temperature of the condensate on the density 
profiles were studied in [34{. The equilibrium properties 
of Newtonian self-gravitating Bose-Einstein condensates 
with short-range interactions were investigated in detail 
in [H, HU ■ The study of the cosmological implications of 
the Bose-Einstein condensation has also become an ac¬ 
tive field of research I n particular, in 143] it was 

shown that condensate dark matter effects can be seen in 
the CMB matter power spectrum if the mass m x of the 
condensate particle lies in the range 15 meV < m x < 700 
meV, leading to a small, but perceptible, excess of power 
at large scales. The possibility of the existence of the 
BEC stars was also considered [40 , while BEC string 
models were investigated in [4§j . 


One of the central issue in the study of Bose-Einstein 
Condensate dark matter halos is the nature, and prop¬ 
erties, of the dark matter particle. One of the most in¬ 
teresting dark matter particle candidate is the axion, a 
(still) hypothetical elementary particle postulated by the 
Peccei-Quinn theory {49j to resolve the strong CP prob¬ 
lem in quantum chromodynamics. It was shown in [50] 
that dark matter axions do form a BEC, as a result either 
of their self interactions, or as a result of their gravita¬ 
tional interactions. It was also proven, from the study of 
the cosmological perturbations, that axion Bose-Einstein 
Condensate dark matter halos differ from standard Cold 
Dark Matter ones on small scales only. Unlike vortices 
in superfluid He 4 and dilute gases, the vortices in the ax¬ 
ion Bose-Einstein Condensate dark matter are attractive 
[5l| . Therefore a large fraction of the vortices in the ax¬ 
ion BEC may form a single big vortex along the rotation 
axis of the galaxy. The resulting enhancement of caustic 
rings could explain the rises in the Milky Way rotation 
curve, usually attributed to caustic rings. The properties 
of the axionic dark matter have been extensively investi¬ 
gated in jH. 


It is the purpose of the present paper to extend the 
previous investigations of the n = 1 polytropic Bose- 
Einstein dark matter halos, by considering in detail and 
in a systematic way the observationally relevant physical 
parameters. As a first step in this study we consider the 
mass and gravitational properties of the dark matter ha¬ 
los in the condensate model. The gravitational potential 
and energy of the halo, its kinetic energy, as well as the 
observationally extremely important logarithmic slopes 
of the density and velocity profiles are obtained. One of 
the most important method which could discriminate be¬ 
tween the condensate model and other dark matter mod¬ 
els is gravitational lensing. We obtain and present in de¬ 
tail, in an analytical form, the physical parameters (sur¬ 
face mass density, deflection angle and magnification etc.) 
necessary for the analysis and comparison of the obser¬ 
vational data with the predictions of the model. Finally, 
the formulations of the virial theorem are obtained for 
rotating and non-rotating dark matter halos, the virial 
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perturbation equation is derived, and some general sta¬ 
bility conditions are formulated. As an application of the 
scalar virial theorem we consider in detail the problem of 
the stability of a slowly rotating and slightly disturbed 
galactic dark matter halo. By introducing a mathemat¬ 
ical description based on the Lagrangian displacements, 
from the scalar virial theorem we obtain the equation 
describing the time evolution of the perturbations. The 
oscillations frequencies and the stability conditions of the 
Bose-Einstein Condensate dark matter galactic halo are 
obtained in the linear approximation. 

The present paper is organized as follows. In Sec¬ 
tion El the basic equations describing the Bose-Einstein 
condensate dark matter halos are presented. In Sec¬ 
tion IIIII the general gravitational properties of the con¬ 
densate halos are investigated. The gravitational lensing 
by dark matter halos is considered in Section IIVI The 
virial theorem for Bose-Einstein Condensate dark mat¬ 
ter halos and their stability properties are presented in 
Section [V] As an application of the scalar virial theorem 
we consider the problem of the stability of a slowly ro¬ 
tating and slightly disturbed galactic dark matter halo 
is considered in Section EO We discuss and conclude our 
results in Section EIIl 


II. BOSE-EINSTEIN CONDENSATE DARK 
MATTER HALOS 


The transition of the dark bosonic component to Bose- 
Einstein Condensed state is assumed to have taken place 
during the cosmological expansion of the Universe at a 
redshift of around z ss 1400 — 1500 [43[. At that mo¬ 
ment the temperature and the matter density of the Uni¬ 
verse did satisfy the condition given by Eq. j2|). Due 
to the expansion of the Universe, its temperature and 
density further decreased. In the following we make the 
fundamental assumption that the dark matter halos are 
composed of a strongly-coupled cold dilute Bose-Einstein 
condensate, at absolute zero temperature. Hence this as¬ 
sumption implies that almost all of the dark matter par¬ 
ticles are in the condensate. Since the dark matter halo 
is dilute and cold, only binary collisions at low energy 
between dark matter particles are relevant. Therefore, 
independently of the details of the two-body potential, it 
follows that the collisions can be characterized by a sin¬ 
gle physical parameter, the s-wave scattering length l a 
[2j, [8(. Therefore, with an excellent approximation, one 
can substitute the interaction potential with an effective 
interaction term of the form Vi (r 1 — r) = UqS (r 1 — r), 
where the coupling constant Uq is related to the scat¬ 
tering length l a through the relation Uq = Airh 2 l a lwiv, 
where m x is the mass of the dark matter particle |2|, [8j . 
In this approach the ground state properties of the dark 
matter are described by the mean-field Gross-Pitaevskii 
(GP) equation 0, gj, 0]. The GP equation for the dark 
matter halos can be derived from the GP energy func¬ 


tional @, n, H, Uli , 


ew = 


h 2 __ ,2 , Uo 


Ivv^T + ^h/’WI 4 - 


2 m x 2 


m x ip* (r) fl ■ (r x p)ip (f) 

1 


dV- 


2 




— Ek + Ei nt + E rot + E, 


grav 


( 3 ) 


where ip (f) is the wave function of the condensate Bi¬ 
ll] . The first term in the energy functional is the quantum 
pressure (kinetic energy Ek), the second is the interac¬ 
tion energy Ei nt , the third term is the rotational energy 
E ro t, with Q the angular velocity, and p the condensate 
momentum, and the fourth term is the gravitational po¬ 
tential energy, respectively. The mass density of the con¬ 
densate is defined as 


p(r) =m x \ip(r)\ 2 , (4) 


and the normalization condition is iV = f \ip (r)| 2 d 3 r, 
where N is the total number of dark matter particles. 
The variational procedure 


5E[ip\—pS J |V’ (r)| 2 d 3 f= 0, (5) 

provides the GP equation as [2, EBl 

- ^ (r) + m x Vg r av (r) ip (r) 

+U 0 | ip (r)\ 2 ip (r) - nL z tp (r) = pip (r), (6) 

where the Lagrangian multiplier p is the chemical po¬ 
tential, and the confining gravitational potential V grav 
satisfies the Poisson equation 

V 2 V grav = 4t xGp. (7) 

The term QL Z = iH(xd v yd x ) is due to the rotation of 
the system about the z axis with an angular frequency 
U. For a rotating Bose-Einstein Condensate dark matter 
halo, its rotational energy may be approximated as [sj 

Erot = 7^ m x J w 2 (H) (x 2 + y 2 ) \ip{f)\ 2 d 3 r, (8) 

where ui is some function of the angular velocity of the 
rotation f2. 

In the physically important case when the number of 
particles in the condensate becomes large enough, the 
quantum pressure term makes a significant contribution 
to the thermodynamical parameters only near the bound¬ 
ary, and it is much smaller than the interaction energy 
term. Thus, for this case the quantum pressure term can 
be neglected (the Thomas-Fermi approximation). When 
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TV — > oo, it can be proven explicitly that the Thomas- 
Fermi approximation becomes exact 0, HHS]. Therefore 
we obtain 

/i - m x Vg r av (f) - \rn x u > 2 ( x 2 + y 2 ) . 

( 9 ) 

In the non-rotating case besides of the conservation of 
the total energy E = Eq + J m x V grav \ip\ d 3 f , where 


rn x 

" m = lk 


E n = 


tl |V'0(fO| 2 + ^ \ip(r)\ 4 


2 m 


d 3 f 


( 10 ) 


where we have denoted 




( 19 ) 


for notational simplicity. 

Note that if the condensate is rotating at an angular 
velocity fl that is greater than the critical one fieri its 
energy is minimized via a creation of vortices. Hence it 
follows that the phase of the condensate order parameter 
changes by 2 tt around a path that includes vortex lines, 


the GP energy functional provides two more conservation 
laws, the continuity equation 


dp 

i +v ^ = 0 ’ 


(ii) 


where the momentum j is given by 
iH 

j = - (ipvr - r wo, (i2) 

and the momentum conservation equation [54 , '55j , 


dji 01 Ijj- OVgray 

dt dxk P dxi 


i, k = 1,2,3, 


( 13 ) 


with 


II ik — 


h 2 
4 m x 


f dtp dip* 
V dxi dx k 


d 2 ip* 
r dxidxk 


+ p 6 ik, 


(14) 


V x « = —V x Vfi (r) = —(f- f,-) , (20) 
777^ m x 


where ft is a unitary vector along a vortex line, fj is the 
radius vector of a vortex line in the plane orthogonal to 
the unit vector n, and S^ is a two-dimensional delta 
function in the corresponding plane sm. 

In order to write the Euler equation Eq. m in a frame 
rotating uniformly with the angular velocity fI we need 
to add to the right-hand side of the Euler equation the 


centrifugal potential 


Q x r 


/2, and the Coriolis accel¬ 


eration 2 u x f2, respectively, where u is the condensate 
velocity in the rotating frame. Therefore the equations 
of motion of an inviscid Bose-Einstein Condensate, ro¬ 
tating at a constant angular velocity fI in the presence 
of a confining gravitational field are given, in Cartesian 
coordinates, by mm 


where c.c denotes complex conjugation, and the quantum 
pressure p is given by 


P = 



2 nH 2 L 


(15) 


By representing the wave function as tp (r, t) = 
p (r,t) exp[iip the velocity of the condensate can 

be defined as v (r,t ) = (h/m x )Vip(f,t). Then Eq. (fill) 
leads to the mass conservation equation 

^+V-(pE)= 0, (16) 

and taking into account the Thomas-Fermi approxima¬ 
tion, Eq. <H3D gives the basic evolution equation of the 
Bose-Einstein Condensate dark matter halos, 

d ,9 dV grav 

~dt pVi + dxk ( pViVk+p ° ik ' = ~ p dx ' 17 

Formally, Eq. m is identical to the Euler equation 
in standard hydrodynamics, however, the distinctive fea¬ 
ture of the Bose-Einstein condensate is that its flow is in 
general irrotational, V x t? = 0 (55]. After applying 
the V 2 operator on both sides of Eq. ©, the Poisson 
equation becomes 


(18) 


duj 


dp 


d 


P dt ^P £ j kl ^ kUl ~ P\ V 9 rav n 


dXn 


dxj 


1 


llxf 


( 2 !) 

where £jki is the completely antisymmetric unit tensor of 
rank three. 


III. MASS AND GRAVITATIONAL 
PROPERTIES OF STATIC CONDENSED DARK 
MATTER HALOS 


For static (non-rotating) condensates to = 0, the gen¬ 
eral solution of Eq. GED, describing the density distribu¬ 
tion p of the static gravitationally bounded single com¬ 
ponent dark matter Bose-Einstein condensate is given by 

uS 


P (r) = Pc 


sin kr 



( 22 ) 


where p c is the central density of the condensate, p c = 
p{ 0). Therefore the Bose-Einstein condensate dark mat¬ 
ter profile is given by the n = 1 polytropic density profile. 

The n = 1 polytropic Bose-Einstein condensate dark 
matter density profile has a well defined boundary radius 
R, at which the dark matter density vanishes, so that 
p(R) = 0. From this condition we obtain R in the form 


X7 2 p + k 2 p = 0, 
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kR = 7r, which fixes the radius of the condensate dark 
matter halo as 


m a 


R k n \l Gml 


(23) 


x 


For r > R the density is smaller than zero, so that the 
n = 1 polytropic density profile cannot be extended to 
infinity, that is, beyond its sharp boundary. For small 
values of r the condensate dark matter profile can be 
written as 


P = Pc 


1 - 


6 R 2 


r 2 + 


120 R 4 


r 4 + ... 


(24) 


The logarithmic slope asE of the Bose-Einstein profile 
is given by [56j 

a DM (r) = = 1 ~ kr cot (hr). (25) 

din r 

Note that Eq. (1251) is not of the power law form 
otDM (r) ~ r ra , where n is a positive number defining the 
steepness of the power law, as are most of the density 
profiles used in the current astrophysical research. For 
r —> 0, cxdm (0) = 0, while at the surface of the dark halo, 
where kR = ir, cxdm diverges, so that limfc r _j. 7r cxdm = 
—oo. For small values of r the logarithmic slope can be 
written as otDAf(r) ~ k 2 r 2 /3 + fc 4 r 4 /45 + O(r) 6 . 

The mass profile m(r ) = 47r f ( ' p(r)r 2 dr of the Bose- 
Einstein condensate galactic halo is 


/ x 47 ip c 

m(r) = -r^-r 
k z 


sin (kr) 


kr 


— cos (kr) 


= -R 2 rp(r)a DM (r). 

7r 

The total mass M of the condensate is 
M(R) = = 1 Pc R\ 

/c* 3 7r 


(26) 


(27) 


which represents a simple cubic proportionality between 
mass and radius. The mass of the condensate is around 
three times smaller than the mass of a constant p c density 
sphere. The central density of the condensate is deter¬ 
mined by the total mass M(R) and the radius of the 
condensate as 


Pc = 


irM(R) 
4 R 3 ' 


(28) 


The mean density (p) = 3M/4nR , 3 of the condensate 
can be obtained as 


/ \ 3p c 3 p c 

{P) ~ K^R? ~ ~7T^~ 


(29) 


Alternatively, we can define the mean density of the dark 
matter halo as 


_ 1 [ K Si(TT) 

P=~5 P\ r ) dr = - 

JL ./n 7T 


p c = 0.5895p c 


(30) 


where Si(i) = f^sin(x)dx/x, and Si(7r) = 1.8519. 

The gravitational potential V grav (r) of the condensed 
dark matter distribution is determined by the equation 


Vgrav(r) = G f 
J r 


m (r 1 ) dr 1 


4 Gp c R 3 . (tt r 

-2- Sln l 

n z r V R 


At small radii the potential behaves as 


r<R. (31) 


_^ r2 + ^ +0 (r)«, (32) 

for r < R. 

The gravitational potential energy U (r) per unit mass 
and inside radius r of the condensed dark matter halo is 
given by 


U(r) = —47 tG j 

Jo 


p(r)m(r) 


2 , 2 GR A p 2 c 

■dr = -— x 


2 nr 


2 + cos 


2 nr \ 

~R~ 


„ ( 27rr 


The total potential energy of the halo is given by 


U(R) = - 


12t r 3 G o 12 Gp\ 


-Pc = 


c 


■R . 


33) 


(34) 


The intrinsic velocity dispersion for an isotropic model 
can be obtained from the definition 


( v r( r )) = 


and is given by 


GM (r') p (r 1 ) dr' 


p(r) 


/ 2 / w 2GR 3 p c . 

W = -5- sm V n 

7 x z r \ R 


/"irr\ 

\-r) 


r /2 


2 GR 2 


(35) 


p(r). (36) 


One can define a kinetic energy K of the halo in terms 
of the average velocity dispersion ay as 

K(r) = | J p(r)ay(r)dV. (37) 


If the velocity dispersion is a constant, we obtain 

sin (kr) 


K ( r ) = ^T<rvr 


and 


k 2 


I<(R) = 


kr 


— cos (kr) 


(38) 


67 T 2 p, 


c 2 


6 p, 


k 3 


_ ~rc J2 p3 

(Jt/ — CT\/kl . 


(39) 


For the ratio of the kinetic and potential energy, we ob¬ 
tain 


K(R) k 2 

|t/(i?OI ~ 2 t tG Pc 


2Gp c R 2 


(40) 
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The tangential velocity of a test particle moving in the 
condensed dark halo can be represented as (Toj 


For the mean of the logarithmic slope of the tangential 
velocity we obtain 


V 2 (r) = 


Gm(r) 47 iGp c 


r 

4 GR 2 
n 


k 2 

p{r)a DM (r). 


sin ( kr ) 
kr 


— cos (kr) 


(41) 


The velocity at the vacuum boundary of the halo has the 
maximum value 


1 P Rcore 

(M = p- / Pv(r)dr 

-ftcore JO 


r-R/1 

_ / r \ 2 

/ 

i-TT 2 (p) 

/ 0 

\RJ 


= K 1 - £)=° ° 88 - 


V 2 (R) 


4 iiGpc 

k 2 


4 Gp c R 2 

71 


GM(R) 

R 


Therefore 

K(R) = cr 2 y = 2 R 2 
117(1?.) I V 2 (R) GM(R) v ' 


(42) 


(43) 


Similarly to the case of the density we can de¬ 
fine a mean value of the tangential velocity as V 2 = 
(1/1?) V 2 (r)dr, and which is given by 


The density corresponding to the inner radius of the 
condensed dark matter halo is 

2 

P (Rcore) — Pc- (50) 

71 

Since the radius R (and consequently also the inner core 
boundary R core ) of the dark matter halo is a universal 
constant, depending only on the fundamental physical 
constants and the mass and scattering length of the dark 
matter particles, it follows that for a given central density 
the product p c R cor e is a universal constant, 


AC 

V 2 = —Si(7i)p c R 2 . 

7 T z 


(44) 


PcR 


core 


7r M(R) 

— ——— = constant. 
8 R z 


(51) 


The logarithmic slope of the tangential velocity is defined 
as Pv = — din V/dlnr, and is given by 


IV. LENSING PROPERTIES OF CONDENSED 
DARK MATTER HALOS 


1+1? Tir cot(7ir/R) — R 


r 2 

7 T / 

r \ 2 

aoM(r) V 

r)J 



(45) 


When r = 0, /3y = — 1, while fdy = 4/2 at the vacuum 
boundary of the dark matter halo. 

One can define the core (inner) radius R CO re of the 
Bose-Einstein condensate dark matter halos as the ra¬ 
dius which satisfies the condition aoM (Rcore) = 1 p5j |. 
By taking into account the definition of the logarithmic 
density slope, we obtain for the core radius the relation 

kRcore = 7t/ 2, Or 


7? 

-Li'C.ore — ~ • 


(46) 


The mean value of the logarithmic density slope within 
the radius 0 < r < R core can be defined as [56j 


(a dm) = 


1 


rRc. 


R, 


OLDM(r) 


core Jo 


= 1- 


2 r /2 


x cot xdx , 


(47) 


and it has the universal value 

(c(dm) = 0.3068. 


(48) 


One of the important ways one could test the Bose- 
Einstein condensate galactic dark matter model is by 
studying the light deflection by galaxies, and in particular 
by studying the deflection of photons passing through the 
region where the rotation curves are flat. In the present 
Section, we consider the basic lensing properties of the 
condensate dark matter halos. For a discussion of the 
gravitational lens properties of scalar field dark matter 
haloes see [57j . 


A. Surface mass density 


The surface mass density of a spherically symmet¬ 
ric lens is obtained by integrating along the line of 
sight of the three-dimensional density profile, E (£) = 
f_°° p (£, r) z, where £ is the radius measured from the 

center of the lens and r = i/£ 2 + z 2 . The surface mass 
density can be written as an Abel transform [ll| , so that 


/ +oo 


p (r) rdr 

/ Vr 2 -£ 2 


(52) 


By inserting the density profile of the Bose-Einstein con¬ 
densate dark matter halos we find 


d rR 

E (£) = 2p c — / 

* 7 f 


sin (kr) dr 
? sjr 2 - £ 2 


(53) 
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Now, introducing a new variable i9 by means of the trans- In this approximation the central surface mass density 
formation r = £cosh$, we obtain can be evaluated as 


jd /*arccosh (R/%) 

E(£)=2p c — / sin [fc£ coshi9] M. (54) 

n Jo 


By taking into account the expansion 
“ r-n s o-(2s+i) 

Sin (fc£ COSh 0) = £ L_Lr_- k 2s+1^2s+l 


s =0 


2s+l 

x E C? s+1 e( 2s - 2l+ W, (55) 
1=0 


where Cf = s\/l\ (s — /)!, we obtain 


S (£) _ R 


°° 2S+ 1 , \S Q-(2s+l)/^<2s+l 

Y' V Z _ L l j 2s+l<r2s+l 

2 Pc n ZJo fco ( 2s + 1 - )! ( 2s “ 21 + 1 '* 


2s-2Z+l 




R? 


2s-2Z+l 


- 1 


(56) 


In the following, for the sake of clarity and simplicity, in¬ 
stead of showing the general results that can be derived 
from Eq. m, we will present only the results obtained 
for some particular values of s. More exactly, we will 
restrict our study of lensing to the specific case s = 7. 
The expansion of the sin x function for s = 7 gives a very 
good approximation of the function for x in the range 
0 < x < 7T. The series expansions can be easily general¬ 
ized and obtained for larger values of s, thus giving the 
possibility of obtaining the lensing parameters of Bose- 
Einstein condensate dark matter halos at any prescribed 
level of precision. In fact, since the condition x < ir 
can be formulated as 7r£cosh0 < 7r, it follows that the 
adopted approximation gives a very good description of 
the condensate dark matter halos for all r < R. 

For s = 7 the surface mass density of the dark matter 
halo is given by 


s(0 

2 p c 


= VR 2 ~ £ 2 



4 

6 \ 

[\ ~ — 

7r 


[\ 18 

+ 600 

35 280 J 


7r 

225 


7r 

~9 

7T 2 

1- 

98 


1 - 

£ 

R 


50 

4 


^ (C 2 
9800/ \R 

El (L 

11025 \R 


or, in numerical form, 


+0.778658 


1.17357- 1.82572 ( ^ 
K 


-0.174402 


(57) 


(58) 


2 p c \ 18 600 35280 ) 


0.5867 R. (59) 


The total mass of the Bose-Einstein condensate model 
can be calculated by integrating the surface mass density 
over the plane of the sky, (R) = 2n / Q K E(£)£d£, and 
is given by 

47T ( 7T 2 7T 4 7T 6 \ , 

~ — f 1-1-) p c R? 

3 V 10 280 15120 ) P 

= 1.2455p c i? 3 , (60) 

which gives an excellent approximation for the exact 
mass formula M(R) = (A/n)p c R 3 = 1.2739 p c R 3 . 

An important quantity for gravitational lensing stud¬ 
ies is the cumulative surface mass density, i.e., the to¬ 
tal mass contained in a infinite cylinder with radius £, 
M (£) = 27t |jE((')f(i(' pl - leof . which for the Bose- 
Einstein condensate dark halo is given by 



Mt(o = p c i 



7T 2 7T 4 

y + 300 



2 (2 7T 2 7T 4 

-7r 2 -1- 

\9 225 14700 

n 2 \ f £\ 4 _ 2tt 6 
11025 ) \RJ 11025 



or, equivalently, 


Af € (f) = 3 . 68689 £ V ^ 2 ~ £ 2 

+0.66349 


1 - 1.5557 | ^ 
K 


-0.14861 


Pc(62) 


B. The lens equation and the deflection of light 

In the following, we consider a simplified gravitational 
lens scenario involving a point source and a circular sym¬ 
metrical lens. The three basic ‘plane’ in thin lens approx¬ 
imation are the source S , the lens L , and the observer O. 
Light rays emitted from the source are detected by the 
lens. For a point-like lens, there will always be (at least) 
two images Si and S 2 of the source. With external shear 
(due to the tidal field of objects outside but near the light 
bundles) there can be more images. The observer sees the 
images in directions corresponding to the tangents to the 
real incoming light paths |58i - l60j . 

In the thin lens approximation, the lens equation for 
axially symmetric lens is [58l - l60} 

V = yE - D LS a , (63) 































where the quantities 77 and £ are the physical positions 
of the source in the source plane and an image in the 
image plane, respectively, a is the deflection angle, and 
Dl, Ds and Dls are the angular distances from observer 
to lens, from observer to source, and from lens to source, 
respectively. By introducing the dimensionless position 
/3 = r//Ds, 9 = £,/Dl and the dimensionless angle a = 
{Dls/D s) d, the thin lens equation can be written as 

P = 0-a(0). (64) 


In the circular-symmetric case the deflection angle is 
given as 


Dlol (£) 


ir § 2 * 

S Jo ^crit 


1 Ms (0 , 

TrYjcrit £ 





(65) 


and can be computed to give 

= Pc 
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714420 
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R 


4?r 4 
11025' 


3 l 5 


7 r 




7T 
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162 
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4tt 6 

893025 


C 2 ^ 


3780 ) 


R J In 


2 R 


R 2 


R + yjR 2 - £ 2 


•.(69) 


In a numerical form we have 


where the convergence s is defined as n (£) = E (£) /E cr ,t, 
where E cr .jt is the critical surface density defined as 
T, crit = c 2 Ds/DkGDlDls- The central convergence, 
k c = k (0) (58H60|, a parameter that determines the lens- 
ing properties of the Bose-Einstein condensate dark mat¬ 
ter halo profiles, is given by the relation 


1.17357p e i? _ 1.17357? tM(R) 
S C rit 4i? 2 S cr ;t 


( 66 ) 


The dimensionless surface mass density can be given by 


4(0 = Pci Vr^^r 2 


- 0.63456 + 0.368622 


-0.159404 


0.0331881 


-0.00430621 


R 


x 1 - 0.62478 In- 


R / 

- 0.634567? 3 x 
2.R \ } 


R+y/R 2 -?, 


R 


(70) 


k{D l 6 ) 


y/R 2 - D'fd 2 
' c 1.17357 R 


1.17357- 1.82572 



C. The magnification factor 


+0.778658 



0.174402 



(67) 


Gravitational lensing effect preserves the surface 
brightness but it causes variations in shape and solid an¬ 
gle of the source. Therefore, the source luminosity is 
amplified by a magnification factor p,, given by [5814601 ] 


For a spherically symmetric lens that is capable of 
forming multiple images of the source, one sufficient con¬ 
dition is k c > 1. In the case k c < 1, only one image 
of the source is formed. Similarly to the case of other 
non-singular profiles, such as the Einasto profile ; 60], the 
Bose-Einstein condensate dark matter density profiles are 
not capable of forming multiple images for any mass. In¬ 
stead, the condition k c > 1 sets a minimum value for lens 
mass required to form multiple images. 

The deflection potential 4 (£) for a spherically sym¬ 
metric lens is given by [58 h60[ 


P = 


1 

(1 — k ) 2 — y 2 


(71) 


where 7 = 7 (£) is the shear. For a spherically symmetric 
lens the shear is given by 


7 (£) = «-« = 


S (£)-£(£) 


E 


crit 


(72) 


where k = E (£) /E cr u, and E (£) is the average surface 
mass density within £ given by 


4(0 = 2 ^ r «(0 


( 68 ) 


= mm 7 - 


(73) 
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By using Eq. m we obtain for E (£) the expression 


S(£) = 2iip c {|) j0.198228 + i/l - x 


0.487671 ( - 0.384069 + 


0.114005 


R 


0.019378£ 8 

7C 


- 0.198228 


(74) 


Therefore for the shear q(£) of a BEC dark matter halo 
we obtain the expression 




- 0.68590- 


0.198228 ( + 0.155028 (j^j - 


0.664 65 


0.198228 


1.4417 


(Z/R) 2 Jl-(t/R ) 2 


(75) 


The magnification by a spherically symmetric lens can 
be written as i58i-l60l 


1 


P = 


(1 — k) (1 + K — 2/v) 


(1 - E(£)/E crit) (1 + E(C)/E crit - 2E(0/Ecrit) ' 

(76) 


The magnification may be divergent for some image 
positions. The loci of the diverging magnification in the 
image plane are called the critical curves. From Eq. d) 
we see that the lensing profile has two critical curves. 
The first curve, 1 — k = 0 is the tangential critical curve, 
which corresponds to an Einstein ring with a given Ein¬ 
stein radius. The second curve, 1 + R — 2k = 0 is the 
radial critical curve, which also defines a ring, and its 
corresponding radius [58j-[60(. 


V. THE VIRIAL THEOREM AND THE 
PERTURBATION EQUATION FOR 
BOSE-EINSTEIN CONDENSATES 

One of the basic relations in theoretical physics, the 
virial theorem, represents a very powerful method for 
the study of the equilibrium and perturbation properties 
of fluids, including the quantum ones. The virial theo¬ 
rem was originally proved, and used for the study of the 
equilibrium and of the stability of rotating fluid bodies, 


bound by self-gravitation (6ll463(. One of the important 
forms of the virial theorem, the tensor-virial theorem, 
reduces the local hydrodynamical Euler equations into 
global virial equations, which give the essential informa¬ 
tion on the structure and stability of the whole gravi¬ 
tating system. An important application of the virial 
theorem is for the study of the small perturbations of in¬ 
compressible uniform ellipsoids when perturbed from an 
initial equilibrium state. In this case, in the absence of 
viscous dissipation, each perturbed virial equation pro¬ 
vides a different set of normal modes [55! • Moreover, in 
the physically very interesting case when the equilibrium 
of the considered system is sustained by an external con¬ 
fining potential, the perturbations of the fluid obviously 
do not modify the confining potential itself. Therefore 
it follows that the tensor virial methods, and the virial 
theorems, can be extended to nonuniform compressible 
flows. The virial theorem represent an extremely power¬ 
ful method for the study of gases with polytropic equa¬ 
tions of state PlU- It has also been extensively used 
in the study of Bose-Einstein Condensates [53l - [55| . 


A. The scalar virial theorem 


In order to derive the scalar virial theorem we consider 
the behavior of the physical parameters of the dark mat¬ 
ter halos under the scaling transformation r —> ar, where 
a is a constant. Then the normalization condition giving 
the total number of particles, N = f \%/j (r)| 2 d 3 r requires 
that ijj (r) —> a~ 3 / 2 ^{r) [1, 0- Thus the total energy 
scales as 


E [a] — a 2 Ek + ot 2 E ro t + a 3 Ei n t + a 1 E grav . (77) 

Since the energy is stationary for any variation of the 
wave function ijj around the exact solution of the Gross- 
Pitaevskii equation, by requiring the energy variation to 
vanish at first order in a, that is, (SE [a|/^a)| a _ 1 = 0, 
we obtain the virial theorem in the form [2j, [§] 

2 Ek — 2 E ro t + 3 Ei n t + E grav = 0. (78) 

By multiplying both sides of Eq. m with p{r) and 
integrating over r we obtain 


pN — E rot + 2 Ei n t + 2 E grav , (79) 


where the interaction energy Ei nt and the gravitational 
energy E grav are defined as E int = (U 0 / 2) J v p 2 (f)d 3 f , 

and V g r a v = (^x/ 2 ) fv V 9 rav (C) P W respectively. 
On the other hand from Eq. © we obtain 

p = —p c + m x V gra v( 0). (80) 

m x 

By using Eqs. d79ll and dSOll we obtain for the total energy 
of the dark matter halo the expression 



U 0 


m^ 


Pc 4” ^x^9 rav (^) 


-Er 


(81) 
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By using the virial theorem we can easily find the con¬ 
dition for the validity of the Thomas-Fermi approxima¬ 
tion. By using the density for the static condensate given 
by Eq. (1221), it is easy to estimate Ek ~ 2 M(R)k 2 /m^ 
and Ei n t ~ UoM(R) 2 k 3 /to 2 , respectively. The Thomas- 
Fermi approximation requires Ek « Ei nt , giving 


N = 


M 


» 


1 

kin 


R 

^Tin 


R » 


to-1 


4 laPc 


(82) 


Hence for systems with enough high particle numbers the 
Thomas-Fermi approximation is always valid. 


B. The tensor virial theorem 


In the following we assume that the dark matter gravi¬ 
tational condensate rotates like a rigid body, and its rota¬ 
tion is supported by an array of vortices. Moreover, it is 
confined by a gravitational field with potential V grav , as¬ 
sumed to satisfy the Poisson equation Eq. ©■ To obtain 
the tensor virial equation for the BEC halo we multiply 
both sides of Eq. m by Xi and we integrates over the 
volume V of the dark matter halo. Thus we first obtain 

[El [Ml 


is the potential tensor of Chandrasekhar [6ll-i63j. By 
defining the kinetic energy tensor as 

Tij = 7}j pUiUjd 3 ?, (89) 

the tensor virial theorem for a rotating Bose-Einstein 
Condensate dark matter halo in a gravitational field can 
be written as [6ll464| 

d f 

.1 P x i u id 3 i' = 2 Tij + SijYl + + Q 2 Rj — 

^jkik^k T 2^ jkl^k I PXiUld V. 

Jv 

(90) 

For a dark matter halo in hydrostatic equilibrium Ut = 
0, and from Eq. m we find 

i’ij + 5ijfi + n 2 Iij — Qjlikflk = 0. (91) 

For i ^ j we have Ty + H 2 /^- — fljlikltk-, while con¬ 
tracting over i and j gives 

<f> + 3n + n 2 I - Qilij Qj = 0, (92) 

where <I> represents the gravitational potential energy. 


d_ 

dt 


[ pXiiijd 3 r = f pxi-d-d 3 r+ f pUiUjd 3 f. (83) 
Jv Jv dt J v 


1. Gravitational energy of a rotating dark matter halo 


The last term on the right-hand of Eq. m can be rewrit¬ 
ten as 


( d 


/ P X '' R~ 

Jy OXj 

f!xr 


d 3 f=n 2 i ij -n j i ik n k , (84) 


where 


Aj — A* — 


pXiXjd 3 r, 


(85) 


is the moment of inertia tensor of the dark matter con¬ 
densate. Assuming that the pressure becomes a constant 
Po outside the dark matter distribution, one can write 


- / x, 


d P , 3 


dx-, 


d 3 r = — Xi.podT,j + / 8 ijpd 3 r = 


$ij / {p-Po)d 3 f=8ijU, (86) 


where Sij is the Kronecker delta symbol. By tak¬ 
ing into account the explicit form of the con¬ 
densate self-gravitational potential, cf> (r, t) = 
— G f v d 3 r'p ( ft ) / | f— 7*1, one can show that 


P x il¥r d 3 r= ~\j v P B v d3 r = -&i . 


dxj 


(87) 


where 


Bij (r, t) = -G f pif ,t) 
Jv 


{Xi - x\) ( Xj - x'j) 


I r — r'l 3 


dV, (88) 


As an application of the tensor virial theorem we de¬ 
termine in the following the gravitational energy of a ro¬ 
tating dark matter halo. If the u-axis is chosen as axis of 
rotation, Eq. (1M1) takes the form 


$ + 3II + n 2 I n = 0, (93) 

where In = f v p (x 2 + y 2 ) d 3 r. For the case of the rotat¬ 
ing Bose-Einstein condensate dark matter halo the equa¬ 
tion of the hydrostatic equilibrium takes the form 


PP 

P 


2V— — V ( V grav + 2 


fl x r 


and can be integrated to give 

1 


2 P — P I Vgrav + 


fI X f 


- 1 / (0) 
v grav 


(94) 


(95) 


where Vgrlv is the gravitational potential at the pole of 
the dark matter halo. Integrating Eq. (K)5l) over the dark 
matter halo volume we obtain 

2n = -2$ + in 2 /o - (96) 

where <F = — (1/2) f pV grav d 3 r. Hence from Eq. (l93l) we 
obtain (65| 

* = \ ~ 3 V g W v M^j . (97) 

For a static halo H = 0, V d P av = GM/R , and we reob¬ 
tain the well-known result $ = —(3/4 )GM 2 /R. 
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C. The perturbation equation 


We consider now the small oscillations about equilib¬ 
rium of a Bose Einstein Condensate dark matter halo. 
We assume that in the stationary state there are no fluid 
motions. In order to obtain the perturbation equation 
we use the tensor virial theorem. Considering only peri¬ 
odic oscillations with frequency w, representing the most 
interesting case from a physical point of view, the La- 
grangian displacement of a mass element dm can be writ¬ 
ten as (r , t) = |*(r) exp ( lot ) [55j, [filT - liul . Due to the 
equation of continuity for such Lagrangian displacements 
dm is a constant, and dp/p + V • £ = 0. By taking the 
Eulerian variation of the tensor virial equation for the 
condensate dark matter halo we obtain 

5— I pXiUjd 3 r = 25Tij + dijSTl + 5$ij + Cl 2 Vij 

dt Jy 

-t- ‘IcjkiQkS I pXiiiid v, (98) 

Jv 

where 

Vij = dlij = I P {Z,iXj + £jXi) d 3 r, (99) 
Jv 


is a tensor symmetric in its indices. By defining its non- 
symmetric part as Vyj = f v p£iXjd 3 r, we have Vij = 
Vi-j + Vj-i. For the variation of we find the equation 



which gives the important first order change in the grav¬ 
itational potential energy, due to a small perturbation of 
the matter in the condensate dark matter system. There¬ 
fore, when there are no condensate dark matter motions 
in the unperturbed frame, the second order virial equa¬ 
tion takes the form Q5. 6lM63l| 

^2 J = ^SiklQl ^ + d<f>ij — 12il2fcVfcj + dijdfi. (101) 

For small periodic perturbations we obtain [55| 

a ~Vi;j = 2£iklQl ( xVk;j + d&ij — 12jl2fcI4y + <5,y<5fl. (102) 

Eq. (11021) contains all the second-harmonic modes of 
the rotating condensate dark matter halo confined by its 
own gravitational field jSjJ iMl - iffl . 


VI. APPLICATIONS OF THE VIRIAL 
THEOREM TO BEC DARK MATTER HALOS 

As a simple example of the application of the virial the¬ 
orem for the study of the stability of the Bose-Einstein 


Condensate dark matter structures we consider the case 
of the slowly rotating and slightly distorted galactic ha¬ 
los. By employing the Lagrangian variables and the 
scalar virial theorem, and by using the linear approxi¬ 
mation, we obtain the angular frequency a of the lowest 
radial mode, as well as the condition for dynamical insta¬ 
bility, depending on the numerical value of the adiabatic 
exponent 7 = 1 + 1 /n, where n is the polytropic index. 


A. The scalar virial theorem for axisymmetric 
BEC dark matter halos 


We consider the BEC dark matter halo as a uniformly 
rotating, homogeneous and compressible spheroid. The 
galactic halo rotates at t = 0 with an angular velocity 
fl (to) = Ho- We restrict our analysis to axisymmetric 
pulsations, and work in a Cartesian coordinate system 
(xi, X 2 , £ 3 ), which rotates at any instant t with the in¬ 
stantaneous angular velocity 12(f) of the galactic halo. 
Moreover, we assume that 12 is oriented along the £3 
axis. We will also ignore dissipative and electromagnetic 
effects. Under such assumptions, the equation of motion 
of the galactic halo, Eq. m, can be written as [6ll - l66| 


where i = 1,2,3 and we have taken into account that 
the angular velocity 12 explicitly depends on time. For 
axisymmetric motions all terms in Eq. (11031) describe mo¬ 
tions in meridional planes, except the_last two terms on 
the left-hand side of the equation j64|, [65}. Therefore 
Eq. (11031) can be split as 


dui 

dt 


1 dp 
p dxi 


dV m 


dxi 


+ (1 - s i3 ) n 2 xi, 


(104) 


where i = 1, 2, 3 and 


tik 3 ( 212Wfc T ^ Xk J 0, 


(105) 


respectively. Equation (11051) can be immediately inte¬ 
grated to give 


[12 {x\ + x 2 )\ = 0. (106) 

Equation (11061) is the equation of conservation of the 
angular momentum for each fluid element, a property 
which is true for axisymmetric motions only. Now we 
multiply Eq. (1 1041) by Xi, we integrate over the total mass 
M of the galactic halo, and add the resulting equations. 
Hence we obtain the scalar virial theorem for the BEC 
halo in the form 0, [65} 


2 dt 2 


f , [ dx k dx k , 

Im Jm dt dt 


<f> 


+3 f —dm + j 12 ( x 2 + £ 2 ) dm, (107) 
Jm P Jm 
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where 


1 f f dmdm' 

2 JmJmV^V 


(108) 


is the gravitational potential energy of the dark matter 
halo. At time t = 0 the dark matter halo is in relative 
equilibrium, and the virial theorem Eq. (11071) gives 


$0 + 3 / —dm + f Qq ( a\ + a 2 ) dm = 0, (109) 

Jm Po Jm 

where the subscript zero indicates the initial values, and 
ai and a 2 represents the t = 0 values of the coordinates 
X\ and X 2 , respectively. 


B. The equations of motion in Lagrangian 
coordinates 


It is more convenient to study the axisymmetric motion 
of the galactic dark matter halos in Lagrangian coordi¬ 
nates, in which all variables characterizing the BEC sys¬ 
tem, Xi, /?, p, and V are expressed as functions of the in¬ 
dependent quantities m and t [65j, [(36}. Therefore the Eu- 
lerian coordinates aq can be expressed as Xi = x. t ( m,t ). 
Hence, when axisymmetry is preserved, there are expand¬ 
ing and contracting solutions for which the coordinates 
of a particle are functions of their initial values and of 
the time only. 

In Lagrangian coordinates the equations of motion of 
the BEC dark matter halo Eqs. (11041) can be reformulated 
as [66| 


dx i d 2 Xi 
dm dt 2 


1 dp 
p dm 


dV c 


grav 


dm 


+ (1 - S i3 ) Vl 2 Xi 


dxj 

dm ’ 


( 110 ) 


where * = 1,2, 3. 

The conservation of mass imposes the condition 


pJ = po, (HI) 

where J is the Jacobian of the coordinate transformation 



J = 


dx\ dx 2 dx 3 
dm dci 2 da 3 


( 112 ) 


As for the pressure of dark matter halo, we will assume 
that it is given by Eq. m, that is, by a poly tropic equa¬ 
tion of state with polytropic index n = 1. We assume 
that the pressure vanishes on the moving surface, at the 
vacuum boundary of the dark matter halo. Moreover, 
the gravitational potential must be continuous across the 
galactic boundary. Once these conditions are satisfied, 
axisymmetric motions are then entirely determined once 
an initial velocity distribution is prescribed at every point 


C. Pseudo-radial oscillations of BEC dark matter 
halos 

In the following we restrict our analysis to the case of 
pseudo-radial pulsations of slowly rotating and slightly 
distorted galactic poly tropic dark matter halos. For 
slowly rotating spheroids in the first approximation one 
can adopt for the evolution of the coordinates a linear 
expression, so that 0 


Xi = m((t), *=1,2,3. (113) 


The equation governing the time evolution of w{t) can 
be obtained from the scalar virial theorem Eq. (11071) in 
Lagrange coordinates. First of all, by using the adopted 
representations of Xi we obtain immediately the identity 


2 dt 2 


f , f dx k dx k .d 2 C 

/ XkX k dm- —rr—rr drn = ^-rP2 I o, 

! m J m dt dt dt 


(114) 


where Iq denotes the initial moment of inertia of the ro¬ 
tating galaxy, and we have used the constancy of the 
mass element dm as we follow the motion of the parti¬ 
cles. In order to obtain the time variation of the gravita¬ 
tional potential energy, given by Eq. (11081) . we take into 
account that \r — r'\ = £( t ) |a — a'|, in which a = r(0) 
and a' = 7^(0) [66]. Then we obtain for the gravitational 
potential energy the expression 


C 2 ' 


(115) 


With the use of Eq. (11111) and of the equation of state of 
the BEC dark matter we obtain 


P = J_Po 
P C 3 Po' 

By integrating the above equation, we obtain 


(116) 


3 f —dm = i (|$ 0 | — 21 ^ 0 ), (117) 

Jm P C 3 


where Kq represents the initial rotational kinetic energy, 
and in order to eliminate the integral over pressure in 
Eq. (11171) we have used the equilibrium condition given 
by Eq. (11091) [66j . The conservation of the angular mo¬ 
mentum immediately gives the equation 

= (US) 

and, consequently, 

[ fl 2 (cc 2 + x 2 ) dm = (119) 

Jm C 


Therefore, with the use of the above results the scalar 
virial theorem Eq. (I107D we obtain the time evolution of 
( as [66| 


ft 

dt 2 



Jo C 3 



2 Kq 

Iq 


( 120 ) 
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Eq. (11201) must be solved with the appropriate bound¬ 
ary conditions. It is also important to mention that the 
quantities $ 0 , Kq and Iq refer to the rotating galactic 
halo. Moreover, in obtaining Eq. (11201) we have assumed 
that the halo rotates as a solid body. 


D. Oscillation frequency of BEC dark matter halos 
in the linear approximation 


where we have taken into account Eq. (l27l) giving the 
mass-radius relation of the static BEC dark matter halo. 
For the period of the oscillations we obtain 


T = 


2tt 


7r 


3/2 _ 


1 


15 y/G Pc 

Pc 

10 -24 g/cm s 


= 1.5745 x 10 16 x 


s. 


(128) 


In the following we consider the small amplitude 
pseudo-radial oscillations of the n = 1 polytropic dark 
matter halos. We can see from Eq. (11201) that the point 
£ = 1 defines an equilibrium state of the galaxy. In order 
to to describe linear and quasi-linear oscillations about 
this equilibrium point, we Taylor expand £ near £ = 1. 
In the first approximation we obtain 

£(i) = l + e(i), (121) 


In the case of perturbed BEC dark matter halos with 
an initial rotation, the stability condition reduces to 

|$o| > K 0 , (129) 

or, equivalently, with the use of Eq. (El), 

> 2/nX (130) 


where e(t) <C 1. In this approximation Eq. (11201) becomes 

m 


dt 2 


+ a 2 e = Xe 2 + 0(e 3 ) 


where 


a 2 = 2 - 


l < f>ol 2IW 


h h ’ 


and 



( 122 ) 


(123) 


(124) 


where flo is the initial angular velocity, and Jq 0 the initial 
moment of inertia of the galactic halo. The period of 
oscillations of a slowly rotating slightly disturbed BEC 
dark matter halo is found as 


v^Tr-y/ In 0 ~ V2n^/In 0 
V\3>o\-K 0 ~ Vl*o 1- K 0 
V32n\JIn 0 

y/2^/2 ngJ no + GM 2 /R 


(131) 


By using again the mass-radius relation for BEC dark 
matter halos given by Eq. m, we obtain 


If we neglect the second and third order powers of e, the 
solution of Eq. (I122D is obtained as 


T 


/~^~ rr 3/2 1 _ \ _ 

V 15 \JGp c \Jl + 15Gp c flg/7r 


(!32) 


e(t) = eo sin (at). (125) 


In terms of the Eulerian variables we obtain £, = 
eo Xi sin (uit). 

Equation (11231) gives the stability condition for BEC 
dark matter halos in the linear approximation. If the 
halo is initially non-rotating, we have Kq = 0, and the 
stability condition reduces to 

cr 2 = 2-^ > 0, (126) 

Iq 


a condition which is always satisfied by BEC dark matter 
halos described by a polytropic equation of state with 
n = 1 and 7 = 2, respectively. With the use of Eq. J97| 
we can estimate the initial gravitational energy of the 
dark matter halo as |<f>o| = (3/4 )GM 2 /R, while for the 
moment of inertia we use the expression Iq = 2Mi? 2 /5 
[65j . Then the oscillation frequency of the halo is given 
by 


_ 15GM _ 15G ^ 
“ 8 R 3 ~ ~ 2 ^ Pc 


( 127 ) 


This equation gives the corrections to the oscillations pe¬ 
riod of the halo due to the presence of an initial rotation 


VII. DISCUSSIONS AND FINAL REMARKS 


The n = 1 polytropic Bose-Einstein condensate dark 
matter model is the simplest existing dark matter model. 
All the properties of the dark matter distributions are 
determined by two parameters only: the mass and the 
radius of the dark matter halo. The central density of 
the dark matter is determined uniquely by M(R ) and R. 
For halos not contaminated with baryonic matter both 
quantities must have the same universal value. However, 
the presence of the baryonic matter may increase the size 
of the galactic halo, thus leading to some variations in the 
total mass and radius of the galactic structures. 

Since all the properties of the Bose-Einstein conden¬ 
sates are determined by two observational parameters 
only, once these parameters are known, the physical prop¬ 
erties of the galactic halos can be predicted by the model. 
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IC 2574 



FIG. 1: Comparison of the predicted rotation curve of the 
dwarf galaxy IC 2574 (solid curve) with the observational data 
presented in [56[]. The assumed mass of the galaxy is M = 
1.64 x 1O 1O M0 (the observationally determined value M = 
1.462 x 1O 1O M0), and the radius of the galaxy is R = 12.6 
kpc. 


M81dwB 



FIG. 2: Comparison of the rotation curve of the dwarf galaxy 
M81 dwB predicted by the condensate dark matter model 
(solid curve) with the observational data presented in [56] . 
The assumed mass of the galaxy is M = 0.33 x 10 9 Mq (the 
observationally determined value M = 0.3 x 10 9 Mq), and the 
radius of the galaxy is R = 1 kpc. 


Thus, for example, the tangential velocity of the test par¬ 
ticles in circular orbits around galaxies is determined by 
the universal equation 


V\r) 


GM |"sin(7rr/i?) nR 

~R r [ nr/R C ° S ~R 


(133) 


Once the mass of the galaxy and its radius is known, 
the rotation curves can be obtained immediately. For 
example, in the case of the dwarf galaxy IC 2574, by 
assuming for the condensate component a mass of M = 
1.64 x 10 10 Mq and a radius of R = 12.6 kpc, the rotation 
curve is given by 


written as 


PcR 


core 


0.00827 


( M(R) \ 
\10 10 Mq) 


R 

10 kpc 


-2 

g/cnr 2 , 


(135) 

and for pure condensate dark matter halos this product 
must be a universal constant. However, the condensate 
dark matter model predicts a dependence on both mass 
and radius of the constant. The presence of baryonic 
matter may determine some modifications and variations 
in its numerical value (for a comparison of the theoretical 
predictions with the observational data see M)- The 
condensate dark matter model also predicts a relation 
between p c Rcore and the central convergence n c of the 
halo of the form 


V 2 (r) = 5.6374 x 10 3 


sm(irr/R) 
irr/R 


— cos ■ 


irR 

“r" 


km 2 /s 2 . 


( 134 ) 

The comparison of the rotation curve predicted by the 
Bose-Einstein condensate dark matter model for the 
galaxy IC 2574 with the observational data obtained in 
[561] is represented in Fig. [TJ The comparison of the pre¬ 
dicted rotation curve for the dwarf galaxy M81 dwB with 
the observational data of [56§ is represented in Fig. [2] 


The condensate dark matter model predicts a mean 
logarithmic density slope of the inner core {a dm) = 
0.3068, while the observed value obtained from the study 
of brown galaxies is {a dm) = — 0.29±0.07 [56j] (the minus 
sign appears due to the opposite sign convention adopted 
in the present paper). In [67[ it was shown, by using the 
mass decomposition of rotation curves using cored haloes, 
that the product of the central density p c and the core 
radius R CO re is a universal constant, independent of the 
galaxy mass. In the case of condensate dark matter ha¬ 
los the product of the central density and core can be 


PcRcore — S72 (136) 

Hence the condensate dark matter model is the only 
existing model that makes easily testable predictions , 
without any need for fitting, and hence it can be easily 
falsified by comparing the theoretical predictions with 
the observations. However, despite its remarkable suc¬ 
cesses at the galactic scale, the Bose-Einstein condensate 
dark matter model needs further testing. One of the 
best possibilities in confirming/ruling out the model is 
through gravitational lensing. In the present paper we 
have obtained for the first time all the relevant quanti¬ 
ties necessary for an in depth comparison of the theo¬ 
retical predictions with the observational data in an an¬ 
alytical, easy to handle form. These general analytical 
formulae enable arbitrary precision calculation, as well 
as the study of the asymptotic behavior near the vac¬ 
uum boundary of the condensate. These formulae can 
be used in strong and weak lensing studies of galaxies 
and clusters of galaxies, where dark matter is assumed 
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to be the dominant component. Moreover, a compari¬ 
son between lensing and rotation curve predictions can 
be done easily. 

The virial theorem is a very useful tool to investigate 
the general properties of the astrophysical systems. The 
scalar virial theorem gives a very powerful constraint for 
the validity of the Thomas-Fermi approximation, which 
can be formulated as 


R > 1.581 x 10 3 


TO 


10- 37 g 
- 1/2 


1/2 


la 


10 


-20 


cm 


- 1/2 


10 


-24 


g/cnr 


cm. 


(137) 


The above constrain shows that the Thomas-Fermi ap¬ 
proximation gives an excellent description of the proper¬ 
ties of condensate dark matter halos. By using the tensor 
virial theorem we have derived the perturbation equation 
of the dark matter halos, which can be efficiently used to 
study the stability of dark matter halos under small per¬ 
turbations. 

Finally, we consider the problem of the mass of the 
dark matter particle. The use of Eq. (l23l) allows us to 
make a first estimate of the physical properties of the 
dark matter particle. As a function of the mass and 
scattering length of the particle the radius R of the con¬ 
densate dark matter halo is given by R = t: h 2 l a /Gm x , 

the total mass of the condensate dark matter halo M(R) 
can be obtained as 

/ fc2 1 \ 3/2 4 

M(if) = to (cS|) < 138 > 


while the mean value ( p) of the condensate density is 
given by the expression (p) = 3p c /7r 2 |]j|. Therefore the 
dark matter particle mass in the condensate is given by 

m 


m x 


(n 2 h 2 l a \ 1/3 
^ GR 2 ) 

6.73 x 10~ 2 [l a (fm)] 1/3 [R (kpc)p 2/3 eXl39) 


For l a fa 1 fm and R fa 10 kpc, the typical mass of the 
condensate particle is of the order of m x fa 14 meV. For 
l a fa 10 6 fm, corresponding to the values of l a observed 
in terrestrial laboratory experiments, m x fa 1.44 eV. 

An important method of observationally obtaining the 
properties of dark matter is the study of the collisions be¬ 
tween clusters of galaxies, like the bullet cluster (IE 0657- 
56) and the baby bullet (MACSJ0025-12). From these 
studies one can obtain constraints on the physical prop¬ 
erties of dark matter, such as its interaction cross-section 
with baryonic matter, and the dark matter-dark matter 
self-interaction cross section. If the ratio a m = cr/m x 
of the self-interaction cross section a = 47 d 2 and of the 
dark matter particle mass m x is known from observa¬ 
tions, with the use of Eq. (11391) the mass of the dark 
matter particle in the Bose-Einstein condensate can be 
obtained as [HJ 


m x = 



By comparing results from X-ray, strong lensing, weak 
lensing, and optical observations with numerical simula¬ 
tions of the merging galaxy cluster IE 0657-56 (the Bul¬ 
let cluster), an upper limit (68 % confidence) for a m of 
the order of a m < 1.25 cm 2 /g was obtained in [68| . By 
adopting for a m a value of a m = 1.25 cm 2 /g, we obtain 
for the mass of the dark matter particle an upper limit 
of the order 

m x < 3.1933 xlO- 3 ^^) (o&7l) g 

= 0-1791 x (^j meV.(141) 

This mass limit is consistent with the limit obtained from 
cosmological considerations in [69j. By using this value 
of the particle mass we can estimate the scattering length 
l a as 


l a < = 1.7827 x 10~ 6 fm. (142) 

V 47r 

A stronger constraint for a m was proposed in (70], so 
that <r m S (0.00335 cm 2 /g, 0.0559 cm 2 /g), giving a dark 
matter particle mass of the order 


/ \ -4/5 

m x fa (9.516 x 10“ 38 - 1.6 70 x 10" 37 ) (g 

= (0.053 - 0.093) x (^J meV, (143) 

and a scattering length of the order of 

l a fa (5.038 - 27.255) x 10" 8 fm. (144) 

Therefore the galactic radii data and the Bullet Cluster 
constraints predict a condensate dark particle mass of 
the order of m x fa 0.1 meV. Recent results on the self¬ 
interacting dark matter cross section have been presented 
in 0, m By using collisions between galaxy clusters 
as tests of the non-gravitational forces acting on dark 
matter, from the dark matter’s lack of deceleration in 
the bullet cluster collision one obtains a self-interaction 
cross-section of the order of a dm/ m < 1.25 cm 2 /g (68% 
confidence limit) for long-ranged forces a From the 
observation of 72 collisions a self-interaction cross-section 
iTDif/m < 0.47 cm 2 /g (95% CL) was inferred. 
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